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ABSTRACT 

We report small-scale clustering measurements from the PRIsm MUlti-object Survey (PRIMUS) 
spectroscopic redshift survey as a function of color and luminosity. We measure the real-space cross¬ 
correlations between 62,106 primary galaxies with PRIMUS redshifts and a tracer population of 
~ 545,000 photometric galaxies over redshifts from z = 0.2 to z = 1. We separately ht a power-law 
model in redshift and luminosity to each of three independent color-selected samples of galaxies. We 
report clustering amplitudes at hducial values of z = 0.5 and L = 1.5L*. The clustering of the 
red galaxies is ~ 3 times as strong as that of the blue galaxies and ^1.5 as strong as that of the 
green galaxies. We also find that the luminosity dependence of the clustering is strongly dependent 
on physical scale, with greater luminosity dependence being found between r = 0.0625Mpc and 
r — 0.25Mpc, compared to the r = 0.5Mpc to r = 2/i“^Mpc range. Moreover, over a 
range of two orders of magnitude in luminosity, a single power-law fit to the luminosity dependence 
is not sufficient to explain the increase in clustering at both the bright and faint ends at the smaller 
scales. We argue that luminosity-dependent clustering at small scales is a necessary component of 
galaxy-halo occupation models for blue, star-forming galaxies as well as for red, quenched galaxies. 
Subject headings: cosmology: observations — galaxies: statistics — galaxies: evolution — galaxies: 

high-redshift — cosmology: large-scale structure of the universe — surveys 


1. INTRODUCTION 

Because baryons are a subdominant component of the 
matter density of the universe in ACDM cosmology, 
galaxies form within dark matter overdensities. Precise 
measurements of galaxy clustering can thus be used to 
probe the underlying dark matter structure, and the de¬ 
pendence of clustering on galaxy properties can be used 
to examine the connection between galaxy formation and 
the large scale structure environment. 

Redshift surveys supply increasingly plentiful data, 
and galaxy clustering measurements continue to offer one 
of the best ways to interpret these surveys in the context 
of the ACDM framework. At low redshifts, for example, 
the 2-degree Field Galaxy Redshift Survey (Colless et al. 
2001), the Sloan Digital Sky Survey (SDSS; York et al. 
2000), and the 6-degree Field Galaxy Survey (Jones et al. 
2009) measured hundreds of thousands of redshifts over 
tens of thousands of square degrees out to z ^ 0.2. At 
intermediate redshifts, the wide-held AGN and Galaxy 
Evolution Survey (AGES; Kochanek et al. 2012) mea¬ 
sured tens of thousands of redshifts over almost eight 
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square degrees, and at higher redshifts, surveys such 
as the DEEP2 Galaxy Redshift Survey (DEEP2; New¬ 
man et al. 2013), zCOSMOS (Lilly et al. 2007), and the 
VIMOS-VLT (Very Large Telescope) Deep Survey (Le 
Fevre et al. 2005) have measured tens of thousands of 
redshifts each over helds up to a few square degrees. 

Previous studies have established that the galaxy clus¬ 
tering signal depends on observational quantities such as 
morphology, luminosity, color, and their physical analogs 
such as stellar mass and star formation rate. For exam¬ 
ple, Davis & Geller (1976) measured steeper autocorrela¬ 
tion functions for elliptical than spiral galaxies. Dressier 
(1980) provided evidence that galaxies with more lumi¬ 
nous spheroidal components preferred higher density re¬ 
gions, and White et al. (1988) showed that, in agreement 
with Cold Dark Matter models, galaxies with higher cir¬ 
cular velocities traced high density environments. More 
recently, observations from the local universe out to red¬ 
shifts of z ~ 1 demonstrate that red, passive galaxies 
form are more highly clustered than blue, star-forming 
galaxies (e.g., Zehavi et al. 2005; Coil et al. 2008; Skibba 
et al. 2009) , and that galaxies with higher luminosity and 
stellar mass are more clustered than those with lower lu¬ 
minosity and mass (e.g., Norberg et al. 2001; Coil et al. 
2006; Meneux et al. 2006; Meneux et al. 2009; Coupon et 
al. 2012; Marulli et al. 2013). In particular, the faint and 
bright ends of the luminosity spectrum of red galaxies 
show increased clustering (e.g., Hogg et al. 2003; Eisen- 
stein et al. 2005; Swanson et al. 2008; Zehavi et al. 2011). 
Using photometric redshift surveys, the dichotomy be¬ 
tween quiescent and star-forming galaxies, and the emer¬ 
gence of their differential clustering, has been studied out 
to beyond redshift of z = 3 (Williams et al. 2009; Hartley 
et al. 2010; Hartley et al. 2013). 

Into this context, the PRIsm Multi-object Survey 
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(PRIMUS; Coil et al. 2011; Cool et al. 2013) measures 
^ 130,000 redshifts from z = 0.2 out to z = 1.2 over 
almost ten square degrees. With a larger survey area 
to reduce cosmic variance and more depth than previ¬ 
ous intermediate redshift surveys, PRIMUS allows for 
the measurement of the evolution of galaxy properties 
over this redshift range and of clustering as a function of 
galaxy properties out to z ^ 1. Only the VIMOS Pub¬ 
lic Extragalactic Redshift Survey (VIPERS; Guzzo et al. 
2014) is comparable in its targeting of a similar num¬ 
ber of redshifts out past z ~ 1 over a wide survey area, 
although to a slightly lower depth. 

The first clustering results from PRIMUS (Skibba et 
al. 2014) used auto-correlations to measure the galaxy 
clustering as a function of luminosity and color over pro¬ 
jected scales of 0.1/i“^Mpc < r < 20/i“^Mpc. This 
work extends those results to smaller scales and exam¬ 
ines the luminosity dependence as a function of color 
and scale. We employ the cross-correlation methodology 
of Eisenstein (2003; hereafter E03) to measure the real- 
space, deprojected clustering of 62,106 PRIMUS galaxies 
with respect to ~ L* tracer galaxies drawn from a parent 
population of ^ 545,000 photometric galaxies from the 
imaging catalogs that overlaps the PRIMUS footprint. 
Eisenstein et al. (2005; hereafter E05) previously applied 
this technique to LRGs in SDSS; we now extend this to 
a wider range of luminosities and colors, to smaller phys¬ 
ical scales, and to higher redshifts. 

By cross-correlating our primary galaxies against a 
sample nearly ten times larger, we avoid the Poisson 
noise inherent in the autocorrelation of small subsam¬ 
ples of galaxies. For example, it is less important, in our 
case, that bins in luminosity be of relatively equal size, 
because the secondary sample defines the environmental 
densities around each population in precisely the same 
manner, and so even bins with relatively small numbers 
of primary galaxies have small error bars. This is impor¬ 
tant because we wish to measure the clustering out to 
z = 1, where we must contend with small total numbers 
of galaxies, while still measuring trends over subsamples 
in color and luminosity. 

Precise measurements of small-scale clustering are cru¬ 
cial to understanding galaxy formation. These cross¬ 
correlation measurements bear directly on the connection 
between dark matter halo properties and galaxy proper¬ 
ties predicted and interpreted with analytic and semi- 
empirical models of galaxy formation (e.g.. Peacock & 
Smith 2000; Berlind & Weinberg 2002; Yang et al. 2003; 
Zehavi et al. 2005; Conroy et al. 2006; Skibba & Sheth 
2009; Behroozi et al. 2010; Moster et al. 2010; Hearin & 
Watson 2013; Masaki et al. 2013), by probing the “one- 
halo” term and the transition to the “two-halo” term. 
By examining jointly the relative clustering of color and 
luminosity selected galaxy samples over these scales, we 
can study how galaxy evolution depends on the local en¬ 
vironment. In this paper (Paper I), we present our initial 
results along with interpretations in the context of halo 
models. In a follow-up paper (Paper II), we will address 
the presence of galactic conformity (see Weinmann et al. 
2006; Kauffmann et al. 2013; Hearin et al. 2014; Hartley 
et al. 2015) in the PRIMUS sample. 

The paper proceeds as follows: Section 2 provides a de¬ 
tailed explanation of the cross-correlation methodology. 
Section 3 overviews the PRIMUS survey and presents the 


sample selection for both the spectroscopic and imaging 
samples. We present our results in Section 4 and inter¬ 
pret them and compare them to the literature in Section 
5. We conclude in Section 6 by highlighting our main 
results and suggesting avenues for future investigation. 
Throughout this paper, we use an = 0.3 flat cos¬ 
mology with h = 1. We use AB magnitudes with dust 
reddening corrections applied (Schlegel et al. 1998). 

2. CROSS-CORRELATION METHODOLOGY 

We measure the real-space galaxy clustering using the 
angular cross-correlation methodology described in EOS. 
This requires both a primary galaxy sample with known 
redshifts and a significantly larger catalog of photometric 
galaxies. In our case, the primary sample uses PRIMUS 
galaxies with spectroscopic redshifts (Section 3.4), and 
the secondary sample uses galaxies from the imaging cat¬ 
alogs that overlap the PRIMUS field. The imaging cat¬ 
alog provides a tracer population with which to measure 
the real-space density around each spectroscopic galaxy. 
We calculate the cross-correlation signal for each primary 
by assuming that all imaging galaxies are at the same 
redshift as the spectroscopic galaxy. This allows us to 
calculate passively evolved, k-corrected absolute magni¬ 
tudes (see Section 3.3) and projected distances for each 
imaging galaxy. We then cross-correlate with the spec¬ 
troscopic PRIMUS primary galaxy only those imaging 
galaxies in a fixed luminosity bin, so as as to have a 
tracer population with a well-understood number den¬ 
sity. Because foreground and background galaxies in the 
secondary sample have no physical correlation with the 
spectroscopic galaxy, only those imaging galaxy that are 
at the same redshift will contribute to the clustering sig¬ 
nal; E05 shows gravitational leasing can be neglected. 

2.1. Computational Details 

By assuming spherical isotropy, angular clustering 
measurements can be inverted into a measure of the real- 
space clustering, ^(r), using an Abel transformation (von 
Zeipel 1908). This deprojection, however, requires com¬ 
puting the derivative of the angular clustering. Instead, 
the E03 method eliminates the need to calculate this 
noisy derivative by integrating the correlation function 
over a spherical three-dimensional window: 

A(o) = . . / 47rr^^(r)lU(r, a) dr (1) 

V (a) Jo 

where the weighted volume V{a) = J^47rr^W(r,a) dr 
and W (r, a) is the smoothing window, such that the cor¬ 
relation function is weighted by IT at a distance r, given 
a scale length a. A (a) is the overdensity of objects from 
the imaging catalog; to wit, a measurement of A (a) = 1 
would mean that the density of imaging objects in the 
window W (r, a) was equal to the mean density of that 
imaging tracer population. 

In this paper, we pick the same smoothing window 
as E05, which relative to a spherical Gaussian window, 
reduces the contribution of spectroscopic-imaging pairs 
at small angular distances to the weighted overdensity, 
A(a), thus reducing the contribution of systematics from 
photometric deblending. 

ir(.,a) = ;^exp(-^) 


( 2 ) 
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Fig. 1. — PRIMUS coverage from the six science fields used in this work. Galaxies are shown in color; red, green, magenta, 
and cyan are used only to distinguish between the spatial jackknife regions used for error estimation, and so repeated colors are 
not representative of any physical difference. Each jackknife region has approximately 1,500 primary spectroscopic PRIMUS 
galaxies. Black points show the secondary imaging sample (sparse sampled by 50% for visual clarity.) Tan points are randoms 
distributed in the masked regions. At most, 100,000 random points are shown, but usually many fewer. The masked border 
regions extends beyond the edges of the figures, but are not shown for clarity. In fact, the masked regions are more densely 
sampled than photometric regions. 


EOS shows that the A(a) statistic can be computed as a 
pairwise summation of the spectroscopic and the imaging 
catalogs, as: 


A(a) 


1 


E 

jG{sp} 


1 

Mzj)y 


E (3) 

kG{img} 


P iG{sp} 


where Ngp is the number of primary spectroscopic ob¬ 
jects, 00 (-z) is the real-space density of objects at red- 
shift Zj, and Rjk is transverse distance, and where the 


weighting function G{R) is defined as: 


G{R) 


1 dF 
RdR’ 



r^W{R) 
\/R? — r'^ 


dr. 


( 5 ) 

( 6 ) 


This method allows us to compute a noisy estimate of 
the overdensity, Aj(a), around each and every spectro¬ 
scopic object. These individual values can then be triv¬ 
ially averaged together over any subset of galaxies within 
the full sample. This permits quick measurements of the 
A(a) statistic with respect to many different dependent 
quantities — such as luminosity, color, and redshift — 
without the computational overhead involved in repeat¬ 
edly calculating the angular clustering of different sub- 
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samples. In fact, without additional assumptions, the 
calculation returns (j)oiz)^{a), which is the overdensity 
multiplied by the number density of the tracer popula¬ 
tion. So in order to recover A (a), we need a model of 
the redshift evolution of the number density, which we 
define in Section 3. 

To correct for excluded regions (bright stars, bad pix¬ 
els, edges, etc.), we create a random galaxy catalog in 
which a dense set of galaxies are distributed within all 
survey gaps and around the borders of the imaging cat¬ 
alog masks. These random points only need be assigned 
an angular position on the sky (as opposed to a redshift, 
color, or luminosity). For each primary galaxy, Aj(a) is 
calculated for both galaxy-galaxy pairs and for galaxy- 
random pairs. These terms are summed together, but 
the galaxy-random sum carries an additional weight to 
account for the difference in projected densities between 
the imaging catalog and the random catalog. This is, 
in effect, a volume completeness correction, upweighting 
Aj (a) to account for any masked volume in the spherical 
window. 

We truncate the pair-wise summation at an inner ra¬ 
dius of 5", below photometric deblending is not suffi¬ 
ciently accurate. We must also truncate the summa¬ 
tion at an outer radius; for this, we choose r = 9a. 
This is large enough that the analytic correction to the 
sum is small, since W (r, a) falls as a Gaussian in r/a; 
meanwhile, it is small enough that we can measure the 
cross-correlation function out to an effective radius of 
2/i“^Mpc even at our lowest redshifts. To account for 
the inner and outer summation limits of our calculation, 
we must include in our computation an analytic correc¬ 
tion term that is a function of F{R). The full derivation 
and implementation are given in E03 and E05. 

Lastly, to correct for completeness, each galaxy 
receives the PRIMUS •primary sample completeness 
weight. These weights include both a priori magnitude- 
and spatial density-dependent sparse sampling, details 
of which can be found in Coil et al. (2011) and Cool et 
al. (2013), as well as an a posteriori spectroscopic suc¬ 
cess rate weight — f tar get and f collision, respectively, as 
described in Equation 1 of Moustakas et al. (2013). 

It is useful to be able to directly compare A(a) with 
^(r). For an assumed power-law form of ^(r) oc r'*', E05 
note that the relationship between the two is given by: 

A(a) = (^2 + e(a). (7) 

In this work, we assume that 7 = —2, so A(a) = C(a)/3 « 
^(1.73a). For ease of comparison with the literature, we 
thus define: 

4(r) = A(r/1.73) ( 8 ) 

We report all scale-dependent measurements as a func¬ 
tion of r. Thus, in order to obtain our results from 
0.0625 h“^Mpc < r < 2.0/i“^Mpc, we choose scale 
factors that range roughly from 0.036/i“^Mpc < a < 
1.16/i“^Mpc. Our statistic, Ca(’')) directly compara¬ 
ble to the usual real-space correlation function, ^(r). 

To calculate statistical errors, we use jackknife resam¬ 
pling of 39 spatially coherent regions, defined such that 
each one has ^ 1,500 primary spectroscopic PRIMUS 
galaxies. Each jackknife region is of order ^ 0.25deg^, 


and is divided in RA and Dec, such that all PRIMUS 
subfields have between three and sixteen subregions, as 
shown in Figure 1. The spatial covariance matrix is com¬ 
puted using the same conventions as Zehavi et al. (2011). 
For overdensities, at scales and rj, the covariance 
between the two values is given as follows: 

^ _ 1 _ _ 

Cov(^Az,CAj) = (9) 

where = 39 in this work, and U . is the mean value of 
all subregions at scale r^. 

3. DATA 

3.1. PRIMUS Motivation 

The PRIsm MUlti-object Survey (PRIMUS; Coil et 
al. 2011; Cool et al. 2013) is a wide-area, spectroscopic, 
faint galaxy survey out to z ^ 1. Using a low-dispersion 
prism instrument, PRIMUS obtained robust redshifts of 
^ 130,000 unique objects to an accuracy of ^^/(l + z) ^ 
0.005 over 9.1 deg^ of sky and to a depth of iAB ~ 23.5. 

To meet these needs, PRIMUS designed a low- 
dispersion prism to be installed on the Magellan I 
(Baade) 6.5m telescope at Las Campanas Observatory. It 
measured redshifts with a resolution of i? = A/AA ~ 40. 
This is substantially higher resolution than most photo¬ 
metric redshifts (A/AA ~ 3 — 5) and still more than a 
factor of two better than the best photometric surveys, 
such as COMBO-17 (Wolf et al. 2004)) or ALHAMBRA 
(Moles et al. 2008) . This reduces the average redshift 
errors, Uzl(^ + z), from ~ 3 — 5% to only 0.5% com¬ 
pared to the usual photometric errors. Moreover, the 
redshift completeness in PRIMUS is not color-dependent 
(Cool et al. 2013). These errors are sufficient to mea¬ 
sure large-scale clustering, and by reducing uncertainty 
in the distance modulus, they also better constrain the 
galaxy luminosities — important if we wish to measure 
luminosity-dependent clustering. 

The IMACS camera on Baade has a effective field of 
view of 0.18 deg^, making it a good choice for a wide-field 
survey. Compared to a photometric redshift survey, the 
decision to use a prism meant that otherwise blank sky 
pixels would instead gather spectral information. The 
prism also allowed for ^2,500 objects to be observed si¬ 
multaneously in each pointing, more than a traditional 
grism or grating. This multiplexing on a large telescope 
made possible both the excellent survey depth and the 
high number of measured redshifts. For a more recent im¬ 
plementation of similar methods, see Kelson et al. (2014). 

3.2. Science Fields and Photometric Catalogs 

The PRIMUS target selection and data reduction 
pipeline are described in detail by Coil et al. (2011) and 
Cool et al. (2013). This works uses six of the science 
fields with PRIMUS redshifts. These are the Chandra 
Deep Field South-SWIRE field (CDES, Giacconi et al. 
2001), the European Large Area ISO Survey-South 1 
field (ELAIS-Sl, Oliver et al. 2000), the DEEP2 02'“’' 
and DEEP2 23^’' fields, the COSMOS field (Scoville et 
al. 2007), and the XMM-Large Scale Structure Survey 
field (XMM-LSS, Pierre et al. 2004). 
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Table 1 

PRIMUS Spectroscopic and Photometric Statistics 


Field 

RA“ 

Dec“ 

A ^ 

■^spec 

^spec 

red 

^ green 

N ^ 

^^blue 

7V.= , 

lack 

A ^ 

^imq 

N!> 

iiriQ 

CDFS 

03:32 

-28:54 

1.95 

12481 

2726 

1258 

8461 

8 

2.14 

74000 

COSMOS 

10:00 

+02:21 

1.03 

7918 

1802 

737 

5379 

5 

1.59 

59000 

DEEP2 02hr 

02:30 

+00:36 

0.58 

4849 

1030 

500 

3319 

3 

0.86 

25000 

DEEP2 23/ir 

23:30 

+00:09 

0.67 

4814 

1261 

374 

3179 

3 

0.96 

25000 

ELAIS SI 

00:36 

-43:30 

0.90 

6676 

1909 

823 

3944 

4 

1.41 

39000 

XMM-LSS 

02:20 

-04:45 

2.88 

25368 

5740 

2254 

17374 

16 

15.72 

322000 

Total 



9.05 

62106 

14504 

5946 

41656 

39 

22.68 

545000 


“ Approximate field centers (J2000); see Coil et al. 2011) for details 

^ Area of primary science fields (in deg^); does not include primary objects removed due to location outside of imaging mask 
^ Number of primary PRIMUS objects used in this work 
Number of primary PRIMUS objects in each color sample 
® Number of jackknife regions used to calculate cosmic variance errors 
^ Exact area of imaging mask (in deg^) 

^ Approximate number of imaging galaxies used in cross-correlation calculation 


For CDFS, we used SWIRE photometry (Lonsdale et 
al. 2003). For ELAIS-Sl, we used combined photometry 
from ESO/WFI and VLT/VIMOS (Berta et al. 2006, 
2008). For XMM-LSS, we used the photometry from the 
CFHT Legacy Survey (Coupon et al. 2009) that has 
been reprocessed and published by Erben et al. (2009). 
For the DEEP2 fields, we use their photometry from 
fields 3 and 4, which are the 23^’" fields and 02^’' fields, 
respectively (Coil et al. 2004). For COSMOS, we used 
the April 2009 data release (Ilbert et al. 2009). Figure 
I shows the sky coverage of the spectroscopic and imag¬ 
ing samples for each science field. Table 1 lists proper¬ 
ties of each field for both the spectroscopic and imaging 
samples, whose selections are described in detail in the 
subsequent section. 

3.3. Imaging Sample Selection 

We select our secondary imaging sample from the par¬ 
ent catalogs in two steps. First, as necessary, we apply 
zero-point corrections, convert from Vega to AB magni¬ 
tudes, apply extinction corrections, and remove stars us¬ 
ing star-galaxy flags as given by each survey. See Coil et 
al. (2011) for further details on these corrections, which 
are identical in our science analysis, except for CDFS 
zeropoints, for which no offset was applied. We remove 
any objects that fall outside of our imaging masks. It is 
imperative that the photometry inside our masks be as 
uniform as possible (e.g., surveyed to the same depth) 
and that it not be affected by cosmic rays, bad pixels, or 
saturation due to nearby bright stars. Thus, we choose 
to be relatively conservative in expanding bright star re¬ 
gions and masking areas in which there is uneven cover¬ 
age (e.g., vignetting at survey edges). Our survey region 
boundaries are defined using mauigle (v2.2, Swanson et 
al. 2008). 

Second, for the entire imaging sample, we precompute 
a grid of k-corrections from the entire range 0.2 < z < 1.0 
in increments of Az = 0.02. This reduces the compu¬ 
tational overhead of selecting the secondary luminosity- 
defined tracer galaxies. We assign to each imaging galaxy 
the absolute magnitude it would have at the redshift of 
the spectroscopic galaxy. Then, we apply the appro¬ 
priate k-correction (kcorrect, version 4.2; Blanton & 
Roweis 2007) and passively evolve the galaxy from this 
redshift to z = 0.1, by linearly interpolating between the 
grid of k-corrections and evolving the imaging galaxy as 


Mg,o.i,k+e = Mg,o,i^k + Qiz - 0.1), with Q = 2.04 mag¬ 
nitudes per unit redshift (Blanton et al. 2003). A pho¬ 
tometric galaxy becomes part of the secondary sample 
for a given primary object if Mg^o.i,fc-i-e falls within an 
empirically defined range. An empirical tracer sample is 
defined as: 

^g,0.1,k+e ~ 0-5 < Mg^o.l,k+e < ^g,0.1,k+e + 1-0 (1^) 

where M* = —19.39 (Blanton et al. 2003). Note that, 
unless otherwise indicated, all absolute magnitudes in the 
paper are k-corrected and passively evolved to z=0.1, but 
henceforth, we will suppress the additional subscripts, 
and write the absolute magnitudes more concisely as Mg. 

Choosing a fixed absolute magnitude range for our 
imaging sample defines a nearly uniform population of 
tracers galaxies. The actual spatial density of the imag¬ 
ing catalog at a particular redshift is unknown because it 
depends on the true evolution of a diverse population of 
galaxies that includes a range of spectral and morpholog¬ 
ical types. If our model evolution differs from the actual 
evolution, then our clustering amplitude measurement 
will evolve anomolously in redshift, even if the under¬ 
lying clustering amplitude of our primary spectroscopic 
sample is constant; hence, we do not attempt to constrain 
redshift evolution in this work. Nonetheless, the relative 
clustering amplitudes of different luminosity and color 
subsamples at a given redshift are still exactly correct, 
because the same secondary tracer population is used to 
measure the cross-correlation around all PRIMUS galax¬ 
ies, and these relative amplitudes (e.g., red vs. blue clus¬ 
tering) can be compared exactly across redshifts. 

While we have attempted to reduce the number density 
evolution of our secondary sample, we are constrained 
by completeness limits in the imaging catalogs at high 
redshift. A lower value of passive evolution, Q, for our 
empirical selection would bring the SDSS results into bet¬ 
ter agreement with, for example, the DEEP2 luminosity 
function (Willmer et al. 2006) at intermediate redshift. 
However, using the higher passive evolution measured by 
SDSS allows us to use the faintest possible selection at 
high redshift, while still using ~ L* tracers at low red¬ 
shift. This implies that we are likely cross-correlating our 
PRIMUS galaxies with a slightly more biased tracer pop¬ 
ulation at higher redshifts than at lower redshifts. We 
account for this offset in our calculation of the real-space 
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number densities of our tracer sample. We integrate the 
luminosity function to obtain the number density for our 
empirically-selected tracer population. We choose the 
theoretical luminosity evolution model so that the pre¬ 
dicted number density at z = 0.5 of a 1.5 magnitude 
bin around M* is roughly equivalent given the luminos¬ 
ity functions of SDSS (Blanton et al. 2003) and DEEP2 
(Wilhner et al. 2006). We used the DEEP2 fits for their 
(z) = 0.9 sample with ‘optimal’ weights. This is still nec¬ 
essarily an approximation because the luminosity func¬ 
tions do not agree on the slope of the faint end. Table 2 
gives the values of (^com in increments of Az = 0.05. 

The theoretical evolution model passive evolution 
model is: 


Mgfi,i^k+e — fbfg,0.1,fc + Q{z — 0.1) (11) 

Q = go(l + <Zi(z-0.1)) (12) 

where go = 2.0 and gi = —0.8. 

3.4. Spectroscopic Sample Selection 

We select our spectroscopic sample from PRIMUS’s 
‘primary’ catalog, whose known selection function al¬ 
lows us to create a statistically complete galaxy sample. 
We select only those galaxies with high-quality redshifts 
{Q > 3) in the range 0.2 < z < 1.0. We have tested our 
analysis on only the highest redshift quality, Qd, selec¬ 
tion, and the difference in results is within our statistical 
errors, so we choose to use the somewhat larger sample. 
From comparisons to higher-resolution spectroscopy, the 
fractional redshift error -\- z) < 0.03 for 92.2% 

galaxies and <0.1 for 98% of galaxies. The dispersion 
among the 92.2% with the highest quality redshifts is 
o'z/(l -I- z) ~ 0.0051 (for details, see Coil et al. 2011). 
As noted in Section 2, each galaxy in the primary cata¬ 
log has an individual completeness weight by which we 
weight all of our clustering measurements. We eliminate 
from our primary sample any galaxy that falls outside 
our imaging mask. This affects at most only a handful 
of objects in all fields except CDFS, where it affects ~ 600 
galaxies. Even then, this is only ^ 1% of our primary 
catalog, and so we err on the side of being conservative 
and avoid galaxies where photometric issues might bias 
our results. Lastly, we exclude any primary galaxies for 
which the angular masks account for more than 50% of 
the projected area on the sky at the largest scale. 

In this paper, we report results in terms of the k- 
corrected and passively evolved absolute magnitudes. 
Mg, of the PRIMUS galaxies. We use the quadratic pas¬ 
sive evolution model given in the previous subsection, 
where go = 2.0 and gi = —0.8. Figure 2 plots absolute 
magnitude versus redshift for our primary sample. 

We divide our primary science sample into three rest- 
frame color bins: red, green, and blue. We do so using a 
linear fit to the red sequence, using Mg and u — g color, 
over the full redshift range of our primary sample, and 
over the magnitude range —21 < Mg^k < —17, where the 
red sequence is well-sampled. We fit using k-corrected 
(but not passively evolved) absolute magnitudes. Specif¬ 
ically, the Mg vs. u — g color space is binned in 0.25 
increments in magnitude and 0.05 in color, and the red 
sequence is defined in each magnitude bin as the peak 
of the color distribution. We fit the red sequence slope. 



Fig. 2. — Contours for 50%, 87%, and 99.7% inclusion in 
absolute Mg versus redshift space for the 62,106 PRIMUS 
galaxies used in this work. Galaxies are k-corrected and pas¬ 
sively evolved to z = 0.1. Top: A histogram of the redshifts 
in Az = 0.1 bins. 


2.5 
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Fig. 3. — Inclusion contours, as in Fig. 2, but in the color- 
magnitude diagram of primary galaxies. Solid lines show color 
cuts between the red, green, and blue samples, as defined by 
Eqns. (12) and (13), for a hypothetical galaxy at z = 0.5, 
but actual color bins are redshift dependent. See Table 1 and 
Section 3.4 for selection details. 



and then shift it in color-space to define the other color 
bins. The shift is chosen so that the color cuts roughly 
match those of Skibba et al. (2014), in which magnitudes 
were k-corrected to z = 0 rather than to z = 0.1, so that 
the results can be more easily compared. We then allow 
for redshift evolution in the color cuts using the values 
found by Aird et al. (2012). We obtain the following 
color cuts for the lower bound of the red sample and the 
upper bound of the blue sample: 
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Table 2 

Modeled Densities of Secondary Sample 


z 

4^com ^ 

z 

^com (■2^) ^ 

z 

4^com (■S) ^ 

0.20 

1.2798 

0.50 

1.0741 

0.80 

0.6306 

0.25 

1.2623 

0.55 

1.0155 

0.85 

0.5424 

0.30 

1.2384 

0.60 

0.9499 

0.90 

0.4545 

0.35 

1.2077 

0.65 

0.8779 

0.95 

0.3693 

0.40 

1.1702 

0.70 

0.7999 

1.00 

0.2649 

0.45 

1.1257 

0.75 

0.7171 




“ The values of the comoving densities 4>com{z), which are calculated 
by matching SDSS (Blanton et al. 2003) and DEEP2 (Willmer et al. 
2006) luminosity functions at z=0.5, are in units of 10~^ Mpc~^ . 
Because our measurements are in physical units, when we calculate 
, we divide out the proper density, = (1 + z)^(pcom- 


{u - g)red = -0.035Mg - 0.065z -f 0.848 (13) 
{u - g)biue = -0.035Mg - 0.065z -f 0.648 (14) 

The green sample is the independent subset from the 
upper blue boundary to the lower red boundary. Figure 3 
plots Mg versus u — g color for our spectroscopic sample, 
with the color cuts for a galaxy at z = 0.5 overlaid. 

Having chosen to use a flux-limited sample to enhance 
the size of our primary sample, we need to account for 
the redshift-dependence in out results. We do this in two 
ways, each of which highlights different science results. 
First, we fit the clustering amplitudes to a joint power- 
law model in redshift and luminosity. For each color 
sample, the model takes the form: 


^data 

^A,i 


■model 

■A 



(15) 


where is the best-fit amplitude at the fiducial nor¬ 

malization of z = 0.50 and L = 1.5L* at the ith scale. 
We fit the model individually at each proper scale over a 
dense grid in redshift while allowing only the amplitude 
and Pi to vary, and then we calculate the best-fit using 

2—4 

Xtotai — xh ®^ch grid point in redshift. By this, 

2—0 

we demand that the redshift evolution of the clustering 
amplitude be the same as all physical scales, but we al¬ 
low the luminosity dependence to vary as a function of 
scale. The normalizations at z = 0.50 and L = 1.5L* 
are chosen at convenient locations near the median red¬ 
shift value of the overall primary sample, such that small 
variations in these normalizations have negligible effect 
on the best-fit clustering amplitude We suppress 

the superscript in our plots. 

Second, we present empirical clustering measurements 
without any underlying assumption that the luminos¬ 
ity dependence is a power-law. However, because we 
are using a flux-limited sample, the mean redshifts of 
brighter magnitude bins are higher than those of lower 
magnitude bins by approximately z ~ 0.1 per magni¬ 
tude. Lest we confuse redshift evolution for luminosity 
dependence, we “detrend” the empirical clustering am¬ 
plitudes on a object-by-object basis before binning them 
in magnitude. Specifically, we divide the measured 
by where we use the color-dependent value of a 

calculated using the above model. Thus, the clustering 



Fig. 4.— Modeled overdensity (^^) around red, green, and 
blue galaxies (red circles, green triangles, and blue squares, 
respectively) as a function of scale h~^ Mpc . The secondary 
sample has a redshift-evolving 1.5 magnitude range, as given 
by Eqn. (9). The model is a joint power-law fit to redshift and 
luminosity dependence, where redshift dependence is fixed 
across scales, but luminosity dependence varies independently 
at each scale. The clustering amplitude is normalized at 
z = 0.5 and L — 1.5L*. The top panel plots whereas 
the bottom panel plots in order to better resolve the 

fluctuations with scale. Values on the horizontal axis are off¬ 
set by 7% for visual clarity. All error bars represent jackknife 
resampling over 39 spatially contiguous regions, as shown in 
Fig. 1. 

amplitudes presented in this way once again represent 
the clustering strength of a fiducial galaxy at z = 0.5. 

4. RESULTS 

In Figure 4 (top), we present our best-fit cross¬ 
correlation measurements of the primary PRIMUS 
galaxies with respect to the photometric galaxies in the 
fixed, k-hevolution corrected, absolute magnitude range 
Mg — 0.5 < M < Mg -I- 1. These clustering amplitudes 
correspond to the fiducial value L = 1.5L* and z = 0.5, 
as given in Equation 13. The results cover proper scales 
from r = 0.0625Mpc to r = 2.0Mpc, and we 
split them according to the color samples defined in Sec¬ 
tion 3. The bottom plot presents so as better to 
display small deviations from a power-law. Note that 
Mg « —19.93 at z = 0.5 using our quadratic pas¬ 
sive evolution model. The corresponding best-fit val¬ 
ues for the redshift evolution is ared = —1.90 ± 0.64, 
oigreen = —1.25±0.92, and auue = —2.10±1.15, and the 
corresponding best-fit luminosity dependence are shown 
in Figure 5. We highlight key points in this section, and 
we discuss them in light of theory and other observations 
in the following section. 
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Fig. 5. — The best power-law exponent fit for luminosity 
(/3) as a function of scale for the amplitudes shown in Fig. 4. 
Red, green, and blue galaxies are represented by red circles, 
green triangles, and blue squares, respectively. The functional 
form of the fit is given by Eqn. (14), and the jointly best-fit 
redshift evolution (which is not scale dependent) is ared ~ 
-1.90±0.64, Ogreen = -1.25±0.92, and auu^ = -2.10±1.15. 
Luminosity dependence is strongest at smaller scales for all 
three color samples. 

The real-space cross-correlation function of PRIMUS 
galaxies shows a strong dependence on color. The clus¬ 
tering amplitude of red galaxies is ~ 1.5 to 3 times that 
of blue galaxies and ^ 1 to 1.5 times that of green galax¬ 
ies. Given the statistical errors, this relative cluster¬ 
ing bias between red and blue galaxies and between red 
and green galaxies appears mostly constant with physical 
scale. However, at 250/i“^kpc, where the blue galaxy 
sample displays a noticeable inflection from a power-law, 
there is an increase in the relative bias of red and blue 
galaxies. Table 4 reports the ratios of the red-green and 
red-blue clustering amplitudes as a function of scale. 

In Figure 5, we show the best-fit luminosity power-law 
exponents as a function of proper scale. We find clear 
evidence that galaxies exhibit increasingly strong lumi¬ 
nosity dependence at the smallest scales, regardless of 
color, whereas luminosity dependent clustering is consis¬ 
tent with zero at 250 h~^ kpc and 500 kpc . The red 
and blue galaxy samples again appear to show slight lu¬ 
minosity dependence at 1 h~^ Mpc and above, but the 
statistical errors are large. However, as shown next, a 
power-law model of the luminosity dependence is not suf¬ 
ficient to explain the variation in clustering as a function 
of luminosity. Values for the best-fit and luminosity 
dependence parameter (/3) — along with their respective 
jackknife errors — are given in Table 3. 

Figure 6 presents as a function of magnitude. Lu¬ 
minosity dependence appears more complicated than a 
simple power-law relation, and so unlike in the Figure 4, 
we show the empirical averages of the clustering ampli¬ 
tudes, not fits. However, because we are using a flux- 
limited sample, the mean redshifts of the brighter mag¬ 
nitude bins are higher than the lower magnitude bins by 
approximately z ~ 0.1 per magnitude, and so as noted 
in Section 3, we detrend the galaxies to our fiducial red- 


shift of z = 0.5. Due to Poisson fluctuations and the large 
number of bins created by subdividing in both color and 
luminosity, there are several bins in which the cluster¬ 
ing is measured to be negative. We plot those points at 
10“"^ instead, so that the Icr errors can be viewed on the 
logarithmic scale. 

The results indicate that the luminosity dependence 
is not smooth over the range of luminosities available 
in PRIMUS. While blue galaxies’ luminosity dependence 
appears to be fairly well-approximated by a power-law, 
the red galaxies have had their luminosity dependence 
underestimated by fitting a power-law. In particular, the 
apparent lack of luminosity dependence at 250kpc 
and 500kpc in Figure 5 obscures non-monotonicity 
in the luminosity dependence. On these scales, red 
galaxies show the weakest clustering at roughly L *, with 
strongly increasing clustering at higher luminosities, and 
weaker increases in clustering strength at fainter lumi¬ 
nosities. The green galaxy sample, while noisy, dis¬ 
plays similar clustering strength to the blue sample at 
the smallest scales, but is indistinguishable from the red 
galaxy sample at larger scales. 

5. DISCUSSION 

5.1. Comparison to Past Studies 

Our results present a methodologically unique look at 
galaxy clustering at a fiducial redshift of z = 0.5. In par¬ 
ticular, we use cross-correlations to extend the study of 
luminosity and color dependence down to 62.5kpc, 
which is smaller than other works in the literature that 
look at similar trends (e.g., Hogg et al. 2003, Coil et 
al. 2008, Zehavi et al. 2011, Marulli et al. 2013, Skibba 
et al. 2014). Because our results examine scales that 
are primarily a measure of the one-halo clustering, mea¬ 
surement of the clustering length or clustering slope are 
not readily comparable with those measured using pro¬ 
jected statistics to much larger scales. However, we can 
still qualitatively compare the clustering as a function of 
color, luminosity, and scale. 

In our best-fit amplitude measurements of color- 
dependent clustering in Fig. 4, we find an offset of ^ 
2 to 3 between red and blue samples at L = 1.5L* in the 
cross-correlation amplitude, which would be equivalent 
to an offset of ~ 4 to 9 in the auto-correlation function. 
This is consistent with both low redshift SDSS results 
(Zehavi et al. 2011) and high redshift DEEP2 results 
(Coil et al. 2008). We find that this relative clustering 
offset of red and blue galaxies continues to 62.5 h~^ kpc 
without any noticeable increase in magnitude. In fact, 
the largest relative clustering difference as a function of 
scale occurs at r = 250kpc, which we attribute to 
the fact the transition between one- and two-halo clus¬ 
tering terms occurring at smaller scales for blue than for 
red galaxies, due to their location in less massive halos 
at fixed luminosity. As seen in Fig. 6, this feature ap¬ 
pears driven primarily by the fainter blue galaxies, but is 
still present in all but the highest luminosity bin. Skibba 
et al. (2014) also see a small but significant increase in 
the red to blue bias for their lowest luminosity threshold 
sample at the same scale (see their Fig. 11). 

We also find a monotonic increase in clustering 
strength from blue to red galaxies, such that green galax¬ 
ies have an intermediate clustering strength. In both 
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Table 3 

Clustering and Luminosity Fits by Scale and Color for Fiducial 2 : = 0.5 and L = 1.5L* 



0.0625/1-1 Mpc 

0.125/1-1 Mpc 

0.250/1-1 Mpc 

0.50 h ^ Mpc 

1.00/i 1 Mpc 

2.00/1-1 Mpc 

cmodeL 

Red 

1457 ± 215 

419 ± 35 

107 ± 9 

28.7 ± 3.6 

7.61 ± 1.35 

2.57 ± 0.88 

Green 

973± 211 

313 ± 45 

70.6 ± 13.7 

20.8 ±3.5 

5.78 ± 1.34 

1.96 ±0.87 

Blue 

780 ± 106 

166 ± 23 

27.8 ±6.0 

9.82 ±2.43 

3.77 ± 1.17 

1.52 ±0.93 


Red 

0.46 ±0.16 

0.44 ± 0.11 

0.084 ± 0.075 

0.007 ±0.066 

0.19± 0.13 

0.42 ±0.25 

Green 

0.89 ± 0.22 

0.40 ± 0.16 

-0.12 ±0.16 

-0.07 ± 0.11 

-0.086 ±0.12 

0.096 ±0.14 

Blue 

0.62 ± 0.07 

0.42 ±0.13 

0.03 ±0.16 

-0.01 ± 0.18 

0.18± 0.19 

0.37 ± 0.21 


Fig. 4 and Fig. 6, we see some evidence that the green 
galaxy population may be a transitional population, un¬ 
dergoing quenching via the environment or internally, 
rather than being simply a mixture of overlapping red 
and blue galaxy distributions (e.g., Mendez et al. 2011). 
At all scales, the clustering of green galaxies is statis¬ 
tically indistinguishable, but almost exclusively lower, 
than that from the red galaxies. However, at scales of 
r = 250 h~^ kpc to r = 1 h~^ Mpc, the blue population 
demonstrates a significantly lower clustering signal than 
either the red of the green samples, whereas at the small¬ 
est scales, there is no statistical difference between the 
green and blue clustering signal. Thus, it is possible that 
at large scales, green galaxies are demonstrating cluster¬ 
ing similar to the red galaxies, while at smaller scales 
showing clustering similar to the blue galaxies. A similar 
result was previously found by Coil et al. (2008) at z ~ 1, 
but not by Zehavi et al. (2011) at z ~ 0. Alternately, 
the green population may be substantially identical to 
the red population, but mixed with some blue galaxies. 
Most directly, however, the color-dependent results in 
the projected auto-correlation measured by Skibba et al. 
(2014) in PRIMUS do not appear to show this scale de¬ 
pendence for the green valley (see their Fig. 14). This 
apparent discrepancy could be a result of the combina¬ 
tion of different color bins and different scales. Skibba 
et al. (2014) measure substantial difference between the 
‘redder’ and ‘reddest’ samples, but little difference be¬ 
tween the ‘redder’ and ‘green’ samples. Furthermore, the 
projected statistics mix real-space scales, so the the con¬ 
vergence of blue and green clustering results seen in our 
work occurs at scales not adequately probed by Skibba 
et al. (2014). Thus we cannot reliably discern between 
the two possibilities, not only because we are limited by 
small number statistics (see Table 5) but particularly be¬ 
cause at the smallest scale the overall difference in clus¬ 
tering amplitude between red and blue decreases relative 
to the difference at r = 250kpc. Nonetheless, this 
points to the impact such selection effects can have and 
motivates studying clustering as a function of intrinsic 
physical, rather than observable, properties, such as that 
being done by Mendez et al. (in prep.), which is exam¬ 
ining PRIMUS clustering as a function of specific star 
formation rate and stellar mass. 

As shown in Fig. 5, we find roughly equal luminos¬ 
ity dependence in all of our color samples. For r > 
250 h~^ kpc , Fig. 6 shows luminosity dependence in the 
red population only our brightest bin. The lack of lu¬ 
minosity dependence at scales of roughly r = 250 — 
500kpc confirms other observations at higher red- 
shift. VIPERS (Marulli et al. 2013) also shows a lack of 


either luminosity or stellar mass dependence at this scale 
in their redshift-space correlation function (see their Fig¬ 
ure 3), while finding more dependence at both smaller 
and larger physical separations. This effect is washed 
out in their Wp measurements, so we should not expect 
to see a similar effect in Skibba et al. (2014), and we do 
not. On the other hand, Meneux et al. (2009) report find¬ 
ing no statistically significant luminosity evolution in the 
zCOSMOS survey and instead report a trend with stellar 
mass, particularly at small scales (r^ < 0.3Mpc). 
Nonetheless, their Fig. 10 shows weak (if insignificant) 
luminosity evolution at small scales, a trend our re¬ 
sults confirm. Moreover, their measurement of stellar 
mass versus clustering amplitude in their z = 0.5 — 0.8 
bin (see their Fig. 15) shows a similar ‘pinching’ at 
Tp ~ 0.3Mpc, in which there is no trend in stel¬ 
lar mass at that scale, but there is a trend at both lower 
and higher scales. Moreover, in this redshift bin, Meneux 
et al. (2009) report a one-halo clustering slope that is 
steeper for the higher stellar mass bins, as well as signifi¬ 
cant excess power in the two-halo term. They discuss the 
possibility that these effects are due to known large-scale 
structure in the COSMOS field at z ~ 0.73. If so, then 
our results, which include the COSMOS field, would be 
influenced by the same structure. We tested the effect of 
removing the COSMOS field entirely from our analysis, 
but while there were minor shifts in clustering amplitude, 
the luminosity trend at small scales, for all color samples, 
remained the same. 

The luminosity dependence that we find for red galax¬ 
ies agrees well with Zehavi et al. (2011) at low redshift, 
who find that for large scales, there is little luminosity 
dependence on the red sequence until L > AL*, which 
indeed would enter in our brightest magnitude bin. Like¬ 
wise, at our smallest scales, we see an increase in clus¬ 
tering among red galaxies. But the increase is not as 
substantial as in Zehavi et al. (2011) — the clustering re¬ 
mains lower for the faint red galaxies than for our bright 
red galaxies — and we only find evidence at the smallest 
scales, whereas Zehavi et al. (2011) observe the increase 
in clustering out to « 2 Mpc . On the bright end, 
again, at intermediate redshifts (0.20 < z < 0.44), E05 
found a monotonic scale dependence of the luminosity 
trend across a range of scales from r « 0.22 — 7 h~^ Mpc 
in a sample of LRGs, with smaller scales showing higher 
luminosity dependence. Meanwhile, at scales of r > 
1.75Mpc, the luminosity trend is much less steep 
for galaxies with L < 3L*, while at smaller scales the 
luminosity dependence was similar regardless of the lu¬ 
minosity range — results with which PRIMUS concurs 
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Fig. 6. — Overdensities for the red, green, and blue sample (red circles, green triangles, and blue squares, respectively) 
as a function of luminosity Mg. Bins are each one magnitude in width, and centered on integer magnitudes from Mg = —17 to 
Mg = —21. Significant breaks from a single power-law in luminosity can be seen, particularly for red galaxies. Measurements 
have been normalized to redshift 2 = 0.5; see Sec. 3.4 for details. 


at higher redshift. 

The luminosity dependent results for blue galaxies 
match the higher redshift DEEP2 results of Cooper et 
al. (2006) and Coil et al. (2008) more than the low red¬ 
shift SDSS results, which show little luminosity depen¬ 
dence for blue galaxies (e.g., Hogg et al. 2003). While 
not strictly a clustering study. Cooper et al. (2006), in 
their analysis of local environmental densities using a 
projected third-nearest-neighbor metric, found that blue 
galaxies showed as much of a luminosity trend as red 
galaxies at 2 ^ 1, and Coil et al. (2008) found that 
blue galaxies exhibited luminosity-dependent clustering 
at small scales. Luminosity dependence in DEEP2 red 
galaxies was not conclusive due to larger error bars and 
a smaller effect size, but in at least one bin in projected 
distance {vp ~ 0.15; see their Fig. 9), there is a hint of a 
similar trend. Additionally, Coil et al. (2008) used four 
luminosity thresholds for each color subsample, but the 


brightest and faintest thresholds only differed by ~ 1 in 
median magnitude, and probe the L* regime, where lu¬ 
minosity dependence is weakest. Our results show that 
by exploring a wider range of luminosities, both blue 
and red galaxies show luminosity dependence down to 
2 ^ 0.5. While Fig. 6 appears to support a closer adher¬ 
ence to a power-law dependence on luminosity for blue 
galaxies compared red galaxies, further investigation of 
this is necessary to determine whether it is statistically 
significant. Regardless of the exact form of the depen¬ 
dence, luminosity dependence in the clustering of both 
red and blue galaxy populations should be considered a 
necessary ingredient in modeling galaxy formation. 

5.2. Physics of Galaxy Evolution 

Both our color and luminosity dependent results have 
consequences for understanding the physics of galaxy as¬ 
sembly and evolution. We touch briefly on the most sig- 
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Table 4 

Overdensities Relative to Red Sample from Fits 



0.0625/1“^ Mpc 

0.125 Mpc 

0.250 h"" Mpc 

0.500Mpc 

1.00^ Mpc 

2.00Mpc 

red /^^’blue 
red!green 

1.87 ±0.38 

1.49 ±0.39 

2.53 ±0.40 
1.34 ±0.22 

3.84 ±0.89 
1.52 ±0.32 

2.93 ±0.81 

1.38 ±0.29 

2.02 ±0.72 

1.32 ±0.38 

1.69 ± 1.18 

1.31 ±0.73 


Table 5 

Clustering Results by Magnitude (r^^A) 




Color 

Ngal 

0.0625/1“^ Mpc 

0.125/1-1 Mpc 

0.250/1-1 Mpc 

0.50/1-1 Mpc 

1.00/i 1 Mpc 

2.00/1-1 Mpc 




Red 

716 

4.27 ± 1.84 

6.38 ± 1.18 

9.24 ± 1.50 

9.45 ± 1.76 

5.06 ±2.37 

7.11 ± 3.46 


-16.5 

-17.5 

Green 

372 

1.60 ± 2.55 

2.53 ± 1.88 

8.27 ± 1.93 

8.46 ±2.32 

8.36 ±3.82 

4.34 ± 3.82 




Blue 

4357 

0.63 ±0.46 

1.89 ±0.41 

2.03 ±0.72 

2.11 ± 0.71 

1.90 ± 1.24 

1.64 ±2.77 




Red 

2725 

1.92 ±0.71 

2.82 ±0.72 

5.88 ± 1.19 

6.97 ± 1.31 

5.95 ± 1.49 

8.39 ± 3.36 


-17.5 

-18.5 

Green 

1338 

0.90 ± 1.04 

2.77 ± 1.02 

4.57 ± 1.51 

4.33 ± 1.20 

6.44 ± 1.99 

5.13 ±3.74 




Blue 

12086 

0.86 ±0.40 

0.94 ± 0.53 

1.46 ±0.46 

1.95 ± 0.64 

2.74 ± 1.47 

4.03 ± 3.36 




Red 

5621 

3.89 ±0.82 

4.54 ± 0.63 

5.62 ± 0.69 

7.17 ± 1.20 

6.33 ± 1.34 

8.02 ±3.61 


-18.5 

-19.5 

Green 

2186 

3.16 ± 1.32 

3.69 ± 1.04 

4.12 ± 0.86 

6.70 ± 1.27 

6.93 ± 1.52 

9.01 ± 3.95 




Blue 

15515 

2.11 ± 0.60 

1.85 ±0.58 

1.81 ± 0.62 

3.24 ± 0.80 

2.55 ± 1.52 

7.28 ± 4.00 




Red 

4277 

7.37 ± 1.94 

6.56 ±0.87 

6.75 ± 0.91 

6.68 ± 1.13 

7.32 ± 1.86 

15.07 ± 5.09 


-19.5 

-20.5 

Green 

1640 

3.13 ± 1.29 

4.63 ± 1.20 

4.67 ± 1.12 

4.93 ± 1.38 

5.84 ± 1.79 

11.06 ±3.81 




Blue 

7669 

4.88 ±0.92 

3.18 ±0.74 

2.56 ±0.53 

2.96 ± 0.94 

4.88 ± 1.60 

9.98 ±5.30 




Red 

917 

13.45 ±3.28 

14.36 ±3.91 

11.63 ±2.80 

13.31 ± 3.22 

14.47 ± 3.20 

24.54 ± 7.47 


-20.5 

-21.5 

Green 

312 

10.10 ±4.58 

10.77 ±4.51 

8.50 ±2.43 

5.23 ±2.59 

5.17 ±3.07 

7.36 ± 5.87 




Blue 

1031 

5.56 ± 2.44 

4.88 ± 2.04 

6.03 ±2.57 

6.21 ± 2.19 

7.40 ± 2.81 

20.49 ± 6.96 



nificant here. 

In general agreement with work at both low and 
high redshift, we find a monotonic increase in cluster¬ 
ing strength from the blue, through green, to the red 
galaxy samples (e.g., Hogg et al. 2003; Coil et al. 2008; 
Zehavi et al. 2011). However, as mentioned in the previ¬ 
ous subsection, we see signs that at the smallest scales, 
green galaxies are clustered similarly to blue galaxies, 
while at larger scales, they cluster as red galaxies. This 
could suggest suggest that they are not a distinct popu¬ 
lation, separate from blue and red galaxies, but that they 
are being quenched (Krause et al. 2013). This cluster¬ 
ing significance for the green galaxies at small scales is 
marginal, however, and as noted, the red and blue popu¬ 
lation are less distinct, and so physically, we cannot reli¬ 
ably distinguish it from the alternative for the green val¬ 
ley galaxies, which is that they are caused by scatter from 
the red and blue populations, or that they are dusty star¬ 
forming galaxies. However, the fact that the green galax¬ 
ies have much more similar clustering to red than blue 
galaxies between r = 250/i“^kpc and r = 1Mpc 
suggests that dust is an unlikely to explain most for the 
effect. Specifically, it appears that if the green valley is 
a mix of the red and blue populations, then the major¬ 
ity of the galaxies must be scattering in from the redder 
sample. On the other hand, if quenching is responsible 
for the signal, then satellites may be quenched as they 
enter the virial influence of a larger halo, while central 
galaxies may be only now be reaching a point at which 
mergers or internal processes are shutting off star forma¬ 
tion. Given the parity in large-scale clustering with red 
galaxies, green galaxies are likely to be hosted by similar 
halos; a lower small-scale clustering, although very small 
in this work, agrees with the DEEP2 results of Coil et 
al. (2008), and the clustering signal of red spiral galax¬ 
ies in Skibba et al. (2009) (see their Fig. 10), which was 


smaller than for red galaxies generally at projected scales 
of Tp < 1 h~^ Mpc . 

We also find that luminosity dependence is strongest 
at the smallest scales, and weakest at scales of r = 
0.25 — 0.5 h~^ Mpc, for all our our color samples. How¬ 
ever, the quantitatively similar fits to the luminosity de¬ 
pendence that we report in Fig. 5 may belie different 
physical causes. In the two smallest bins in scale, red 
galaxies begin to show increased clustering (although not 
to the extent of Zehavi et al. 2011) at faint luminosi¬ 
ties, while to statistical significance, the blue galaxies’ 
clustering amplitudes plateau. This should be unsur¬ 
prising if external quenching mechanisms are at work: 
less luminous, blue galaxies will be correlated with star¬ 
forming galaxies with lower stellar mass and less total gas 
mass. They will remain star-forming longer if they live 
in matter underdense environments; if they are satellites 
in more matter overdense regions, physical quenching 
mechanisms will more quickly end star-format ion, mov¬ 
ing them onto the red sequence. Thus, more extant satel¬ 
lites will contribute to the red galaxy clustering signal, 
and the contributions to the blue signal will primarily be 
from the underdense field. While recent galactic confor¬ 
mity observations show that quenched fractions around 
passive hosts are higher than around star-forming hosts, 
these results are not in conflict. Physical quenching 
mechanisms will still move some less-luminous blue satel¬ 
lites onto the red sequence, whether by mass-dependent 
quenching (Phillips et al. 2014) or via feedback (Hartley 
et al. 2015) or some other process. Thus, it is not neces¬ 
sary that all low-luminosity blue satellites are quenched; 
the shape of the clustering signal is expected if lower- 
luminosity blue galaxies quench more easily than more 
luminous blue galaxies. 

Likewise, the relative bias of red and blue galaxies is 
particularly strong at r = 0.25Mpc, where there is 
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an inflection in the clustering of the blue galaxy sam¬ 
ple that is at least ^ 2a stronger than at the other 
scales. These scales roughly correspond to the transi¬ 
tion between the one- and two-halo clustering regimes. 
Zheng et al. (2007) fit HOD models to intermediate red- 
shift DEEP2 observations and find the transition from 
the one- to two-halo term occurs at projected distances 
of Vp ^ 0.4 — 0.6Mpc, with the higher luminosity 
thresholds having larger transition distances. As those 
thresholds range from Mb < —19 to Mb < —20.5 
(and noting that Mg ^ Mb), while the mean luminosi¬ 
ties of the red and blue samples are Mg = —19.10 and 
Mg = —18.67, respectively, we would expect the halo 
transition to occur at a somewhat smaller physical scale. 
Thus, this signal may result from blue galaxies being 
found in less massive halos than their equal-luminosity, 
red counterparts. However, future modeling in mock ha¬ 
los of the cross-correlation statistic is clearly needed to 
confirm this as the likely cause. 

6. CONCLUSIONS 

We have measured the clustering of 62,106 spectro¬ 
scopic galaxies from the PRIMUS survey between z = 0.2 
and 0 = 1.0. We used a method previously used prin¬ 
cipally for bright LRGs (E05; see also Hogg et al. 2003) 
to measure the cross-correlation of roughly L* secondary 
galaxies around PRIMUS primary galaxies. We use this 
method to maximize the statistical power of our primary 
sample, while subdividing it as a function of color, lumi¬ 
nosity, and scale. Our key results are that: 

1. We present precise and detailed measurements of 
strong — but notably scale-dependent — luminos¬ 
ity dependence in the clustering amplitude at scales 
of r = 0.0625 — 2.0 h~^ Mpc . This luminosity de¬ 
pendence is evident separately for each of the red, 
green, and blue galaxy samples, showing that lu¬ 
minosity dependence in clustering is not a color- 
dependent phenomena, nor a side-effect of chang¬ 
ing red fractions in higher luminosity samples. 

2. Luminosity dependence is present over the full lu¬ 
minosity range for r < 0.125Mpc, while at 
larger scales, the luminosity dependence is only evi¬ 
dent at Mg = —21, while the clustering plateaus for 
fainter galaxies. Luminosity dependence is a thus 
a complicated emergent phenomena that traces the 
small-scale effects of galaxy formation and evolu¬ 
tion. Reproducing color and luminosity binned 
clustering, in addition to average clustering, pro¬ 
vides better constraints on the addition of new 
galaxy physics and feedback prescriptions in these 
models. 
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3. We find that red galaxies cluster more strongly 
than green galaxies and much more strongly than 
blue galaxies at an effective redshift of z = 0.5. 
On average, over physical scales of r = 0.0625 — 
2.0 h~^ Mpc, we detect a relative bias of red to blue 
galaxies of ~ 2—3 and a relative bias of red to green 
galaxies of ^ 1.5. This agrees with previous work 
at both lower and higher redshifts. 

4. We detect a maximum relative bias of 3.8 ± 0.9 be¬ 
tween red and blue galaxies at r = 250/i“^kpc. 
We posit that the significantly non-power-law be¬ 
havior of the blue galaxy correlation function at 
this scale that leads to this large relative bias is in¬ 
dicative of the one- to two-halo transition occurring 
at a smaller scales for blue galaxies. 

This paper is the second in a series of papers quantify¬ 
ing the clustering properties of PRIMUS galaxies. It fol¬ 
lows Skibba et al. (2014), which examined the color and 
luminosity dependence of the auto-correlation function, 
Wp{rp), out to Tp = 30/i“^Mpc. Additionally, Mendez 
et al. (in prep.) will measure the auto-correlations as 
a function of stellar mass and specific star formation 
rate. Skibba et al. (2015) analyze the stellar mass de¬ 
pendent clustering of PRIMUS galaxies using analytic 
models and mock galaxy catalogs. Bray et al. (in prep.) 
will use the real-space, cross-correlation statistics used in 
this work to probe three-dimensional galactic conformity. 
This follow-up paper (Paper H) will use both luminosity- 
and color-selected secondary galaxies to examine the red 
fractions around PRIMUS galaxies both in the ‘1-halo’ 
and the ‘2-halo’ regime. 
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